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Abstract 

We have searched for the production of a selectron and a squark in e~^p collisions 
at a center-of-mass energy of 300 GeV using the ZEUS detector at HERA. The 
selectron and squark are sought in the direct decay into the lightest neutralino in 
the framework of supersymmetric extensions to the Standard Model which conserve 
ii-parity. No evidence for the production of supersymmetric particles has been 
found in a data sample corresponding to 46.6 pb~^ of integrated luminosity. We 
express upper limits on the product of the cross section times the decay branching 
ratios as excluded regions in the parameter space of the Minimal Supersymmetric 
Standard Model. 



1 



The ZEUS Collaboration 



J. Breitweg, M. Derrick, D. Krakauer, S. Magill, D. Mikunas, B. Musgrave, J. Repond, 
R. Stanek, R.L. Talaga, R. Yoshida, H. Zhang 
Argonne National Laboratory, Argonne, IL, USA p 

M.C.K. Mattingly 

Andrews University, Berrien Springs, MI, USA 

F. Anselmo, P. Antonioli, G. Bari, M. Basile, L. Bellagamba, D. Boscherini, A. Bruni, 

G. Bruni, G. Cara Romeo, G. Castellini^, L. Cifarelli^, F. Cindolo, A. Gontin, N. Goppola, 
M. Corradi, S. De Pasquale, P. Giusti, G. lacobucci, G. Laurenti, G. Levi, A. Margotti, 
T. Massam, R. Nania, F. Palmonari, A. Pesci, A. Polini, G. Sartorelli, Y. Zamora Garcia^, 
A. Zichichi 

University and INFN Bologna, Bologna, Italy ^ 

C. Amelung, A. Bornheim, I. Brock, K. Coboken, J. Crittenden, R. Deffner, M. Eckert, 
M. Grothe^, H. Hartmann, K. Heinloth, L. Heinz, E. Hilger, H.-P. Jakob, A. Kappes, 
U.F. Katz, R. Kerger, E. Paul, M. Pfeiffer, J. Stamm^ H. Wieber 

Physikalisches Institut der Universitdt Bonn, Bonn, Germany 

D. S. Bailey, S. Campbell-Robson, W.N. Cottingham, B. Foster, R. Hall- Wilton, G.P. Heath, 

H. F. Heath, J.D. McFall, D. Piccioni, D.G. Roff, R.J. Tapper 

H.H. Wills Physics Laboratory, University of Bristol, Bristol, U.K. ° 

M. Gapua, L. lannotti, A. Mastroberardino, M. Schioppa, G. Susinno 
Calabria University, Physics Dept. and INFN, Cosenza, Italy ^ 

J.Y. Kim, J.H. Lee, I.T. Lim, M.Y. Pac^ 
Chonnam National University, Kwangju, Korea ^ 

A. CaldwelF, N. Cartiglia, Z. Jing, W. Liu, B. Mellado, J. A. Parsons, S. Ritz^, S. Sampson, 
F. SciuUi, P.B. Straub, Q. Zhu 

Columbia University, Nevis Labs., Irvington on Hudson, N.Y., USA * 

P. Borzemski, J. Ghwastowski, A. Eskreys, J. Figiel, K. Klimek, M.B. Przybycieh, L. Za- 
wiejski 

Inst, of Nuclear Physics, Cracow, Poland ^ 

L. Adamczyk^, B. Bednarek, M. Bukowy, A.M. Gzermak, K. Jeleii, D. Kisielewska, 
T. Kowalski, 

M. Przybycieii, E. Ruhkowska-Zar§bska, L. Suszycki, J. Zajg-c 

Faculty of Physics and Nuclear Techniques, Academy of Mining and Metallurgy, Cracow, 
Poland ^ 

Z. Dulihski, A. Kotahski 

Jagellonian Univ., Dept. of Physics, Cracow, Poland'' 



I 



G. Abbiendi^°, L.A.T. Bauerdick, U. Bchrcns, H. Bcicr^^, J.K. Bienlein, K. Desler, 
G. Drews, U. Fricke, I. Gialas^^, F. Goebel, P. Gottlicher, R. Graciani, T. Haas, W. Hain, 

G. F. Hartner, D. HaselP^ K. Hebbel, K.F. Johnson^'', M. Kasemann, W. Koch, U. Kotz, 

H. Kowalski, L. Lindemann, B. Lohr, M. Martinez, J. Milewski, M. Milite, T. Monteiro^^, 

D. Notz, A. Pellegrino, F. Pelucchi, K. Piotrzkowski, M. Rohde, J. Roldan^^, J.J. Ryan^"^, 
P.R.B. Saull, A. A. Savin, U. Schneekloth, O. Schwarzer, F. Selonkc, S. Stonjek, B. Surrow^*, 

E. Tassi, D. Westphal^^, G. Wolf, U. WoUmer, C. Youngman, W. Zeuner 
Deutsches Elektronen- Synchrotron DESY, Hamburg, Germany 

B.D. Burow, C. Coldewey, H.J. Grabosch, A. Meyer, S. Schlenstedt 
DESY-IfH Zeuthen, Zeuthen, Germany 

G. Barbagli, E. Gallo, P. Pelfer 
University and INFN, Florence, Italy ^ 

G. Maccarrone, L. Votano 

INFN, Laboratori Nazionali di Frascati, Frascati, Italy ^ 

A. Bamberger, S. Eiscnhardt, P. Markun, H. Raach, T. Trefzger^°, S. Wolfle 
Fakultdt fiir Physik der Universitdt Freiburg i.Br., Freiburg i.Br., Germany 

J.T. Bromley, N.H. Brook, P.J. Bussey, A.T. Doyle^^ S.W. Lee, N. Macdonald, G.J. Mc- 
Cance, D.H. Saxon, 

L.E. Sinclair, I.O. Skillicorn, E. Strickland, R. Waugh 

Dept. of Physics and Astronomy, University of Glasgow, Glasgow, U.K." 

I. Bohnet, N. Gcndner, U. Holm, A. Mcycr-Larsen, H. Salchi, K. Wick 
Hamburg University, I. Institute of Exp. Physics, Hamburg, Germany"^ 

A. Garfagnini, L.K. Gladilin^^, D. Kcira^^, R. Klanncr, E. Lohrmann, G. Poelz, F. Zetsche 

Hamburg University, II. Institute of Exp. Physics, Hamburg, Germany 

T.C. Bacon, I. Butterworth, J.E. Cole, G. Howell, L. Lamberti^^, K.R. Long, D.B. Miller, 

N. Pavel, A. Prinias^^, J.K. Sedgbeer, D. Sideris, R. Walker 

Imperial Gollege London, High Energy Nuclear Physics Group, London, U.K. ° 

U. Mallik, S.M. Wang, J.T. Wu^^ 

University of Iowa, Physics and Astronomy Dept., Iowa City, USA p 
P. Cloth, D. Filges 

Forschungszentrum, Jiilich, Institut fiir Kernphysik, Jiilich, Germany 

J.I. Fleckl^ T. Ishh, M. Kuze, I. Suzuki^^, K. Tokushuku, S. Yamada, K. Yamauchi, 

Y. Yamazaki^^ 

Institute of Particle and Nuclear Studies, KEK, Tsukuba, Japan ^ 

S.J. Hong, S.B. Lee, S.W. Nam^^, S.K. Park 
Korea University, Seoul, Korea ^ 



II 



H. Lim, I.H. Park, D. Son 

Kyungpook National University, Taegu, Korea ^ 

F. Barrciro, J. P. Fernandez, G. Garcia, C. Glasnian"^°, J.M. Hernandez, L. Hervas^^, 

L. Labarga, J. del Peso, J. Puga, J. Terron, J.F. de Troconiz 

Univer. Autonoma Madrid, Depto de Fisica Teorica, Madrid, Spain " 

F. Corriveau, D.S. Hanna, J. Hartmann, L.W. Hung, W.N. Murray, A. Ochs, M. Riveline, 
D.G. Stairs, M. St-Laurent, R. Ullmann 

McGill University, Dept. of Physics, Montreal, Quebec, Canada"', ^ 
T. Tsurugai 

Meiji Gakuin University, Faculty of General Education, Yokohama, Japan 

V. Bashkirov, B.A. Dolgoshein, A. Stifutkin 

Moscow Engineering Physics Institute, Moscow, Russia ^ 

G. L. Bashindzhagyan, P.F. Ermolov, Yu.A. Golubkov, L.A. Khein, N.A. Korotkova, 
LA. Korzhavina, V.A. Kuzmin, O.Yu. Lukina, A.S. Proskuryakov, L.M. Shcheglova^^, 
A.N. Solomin^^ S.A. Zotkin 

Moscow State University, Institute of Nuclear Physics, Moscow, Russia 

C. Bokel, M. Botje, N. Briimmer, J. Engelen, E. Koffeman, P. Kooijman, A. van Sighem, 

H. Tiecke, N. Tuning, W. Verkerke, J. Vossebeld, L. Wiggers, E. de Wolf 
NIKHEF and University of Amsterdam, Amsterdam, Netherlands ' 

D. Acosta^^, B. Bylsma, L.S. Durkin, J. Gilmore, C.M. Ginsburg, G.L. Kim, T.Y. Ling, 
P. Ny lander, T.A. Romanowski^^ 

Ohio State University, Physics Department, Golumbus, Ohio, USA ^ 

H.E. Blaikley, R.J. Gashmore, A.M. Gooper-Sarkar, R.G.E. Devenish, J.K. Edmonds, 
J. GroBe-Knetter^^ N. Harnew, G. Nath, V.A. Noyes^^ A. Quadt, O. Ruske, J.R. Tickner^^ 
R. Walczak, D.S. Waters 

Department of Physics, University of Oxford, Oxford, U.K. ° 

A. Bertolin, R. Brugnera, R. Garlin, F. Dal Gorso, U. Dosselli, S. Limentani, M. Morandin, 

M. Posocco, L. Stanco, R. Stroili, G. Voci 

Dipartimento di Fisica dell' Universitd and INFN, Padova, Italy ^ 

J. Bulmahn, B.Y. Oh, J.R. Okrasihski, W.S. Toothacker, J.J. Whitmore 
Pennsylvania State University, Dept. of Physics, University Park, PA, USA * 

Y. Iga 

Polytechnic University, Sagamihara, Japan ^ 
G. D'Agostini, G. Marini, A. Nigro, M. Raso 

Dipartimento di Fisica, Univ. 'La Sapienza' and INFN, Rome, Italy ^ 
J.G. Hart, N.A. McGubbin, T.P. Shah 

Rutherford Appleton Laboratory, Ghilton, Didcot, Oxon, U.K. ° 



HI 



D. Epperson, C. Heusch, J.T. Rahn, H.F.-W. Sadrozinski, A. Seiden, R. Wichmann, 
D.C. Williams 

University of California, Santa Cruz, CA, USA ^ 

H. Abramowicz^^, G. Briskin^^, S. Dagan^^, S. Kananov^^, A. Levy^^ 

Raymond and Beverly Sackler Faculty of Exact Sciences, School of Physics, Tel-Aviv 

University, 

Tel-Aviv, Israel ^ 

T. Abe, T. Fusayasu, M. Inuzuka, K. Nagano, K. Umemori, T. Yamashita 
Department of Physics, University of Tokyo, Tokyo, Japan ^ 

R. Hamatsu, T. Hirose, K. Homma^°, S. Kitamura"*^, T. Matsushita 
Tokyo Metropolitan University, Dept. of Physics, Tokyo, Japan ^ 

M. Arneodo^^ R. Cirio, M. Costa, M.I. Ferrero, S. Maselli, V. Monaco, C. Peroni, 
M.C. Petrucci, M. Ruspa, R. Sacchi, A. Solano, A. Staiano 

Universitd di Torino, Dipartimento di Fisica Sperimentale and INFN, Torino, Italy ^ 
M. Dardo 

// Faculty of Sciences, Torino University and INFN - Alessandria, Italy ^ 

D.C. Bailey, C.-P. Fagerstroem, R. Galea, K.K. Joo, G.M. Levman, J.F. Martin R.S. Orr, 
S. Polenz, A. Sabetfakhri, D. Simmons 

University of Toronto, Dept. of Physics, Toronto, Ont., Canada"" 

J.M. Butterworth, CD. Catterall, M.E. Hayes, E.A. Heaphy, T.W. Jones, J.B. Lane, 
R.L. Saunders, M.R. Sutton, M. Wing 

University College London, Physics and Astronomy Dept., London, U.K. ° 

J. Ciborowski, G. Grzelak^^, M. Kasprzak, R.J. Nowak, J.M. Pawlak, R. Pawlak, B. Smal- 

ska, 

T. Tymieniecka, A.K. Wroblewski, J. A. Zakrzewski, A.F. Zarnecki 
Warsaw University, Institute of Experimental Physics, Warsaw, Poland ^ 

M. Adamus 

Institute for Nuclear Studies, Warsaw, Poland ^ 

O. Deppe, Y. Eiscnbcrg"''^, D. Hochman, U. Karshon^^ 

Weizmann Institute, Department of Particle Physics, Rehovot, Israel 

W.F. Badgett, D. Chapin, R. Cross, S. Dasu, C. Fondas, R.J. Loveless, S. Mattingly, 
D.D. Reeder, W.H. Smith, A. Vaiciulis, M. Wodarczyk 
University of Wisconsin, Dept. of Physics, Madison, WI, USA ^ 

A. Deshpande, S. Dhawan, V.W. Hughes 

Yale University, Department of Physics, New Haven, CT, USA ^ 

S. Bhadra, W.R. Friskcn, M. Khakzad, W.B. Schmidke 

York University, Dept. of Physics, North York, Ont., Canada^ 



IV 



^ also at IROE Florence, Italy 

^ now at Univ. of Salerno and INFN Napoli, Italy 

^ supported by Worldlab, Lausanne, Switzerland 

^ now at University of California, Santa Cruz, USA 

^ now at C. Plath GmbH. Hamburg 

^ now at Dongshin University, Naju, Korea 

^ also at DESY 

^ Alfred P. Sloan Foundation Fellow 

^ supported by the Polish State Committee for Scientific Research, grant No. 2P03B14912 

now at INFN Bologna 

now at Innosoft, Munich, Germany 

now at Univ. of Crete, Greece, partially supported by DAAD, Bonn - Kz. A/98/16764 
now at Massachusetts Institute of Technology, Cambridge, MA, USA 
visitor from Florida State University 

supported by European Community Program PRAXIS XXI 

now at IFIC, Valencia, Spain 
^'^ now a self-employed consultant 
i« now at CERN 

now at Bayer A.G., Lever kusen, Germany 

now at ATLAS Collaboration, Univ. of Munich 

also at DESY and Alexander von Humboldt Fellow at University of Hamburg 

22 on leave from MSU, supported by the GIF, contract 1-0444-176.07/95 

23 supported by DAAD, Bonn - Kz. A/98/12712 
2^ supported by an EC fellowship 

25 PPARC Post-doctoral fellow 

2^ now at Applied Materials Inc., Santa Clara 

2''' now at Osaka Univ., Osaka, Japan 

2^ supported by JSPS Postdoctoral Fellowships for Research Abroad 
2^ now at Wayne State University, Detroit 

3° supported by an EC fellowship number ERBFMBICT 972523 

3^ partially supported by the Foundation for German-Russian Collaboration DFG-RFBR 
(grant no. 436 RUS 113/248/3 and no. 436 RUS 113/248/2) 

32 now at University of Florida, Gainesville, FL, USA 

33 now at Department of Energy, Washington 

3^ supported by the Feodor Lynen Program of the Alexander von Humboldt foundation 

35 Glasstone Fellow 

3^ now at CSIRO, Lucas Heights, Sydney, Australia 

3''' an Alexander von Humboldt Fellow at University of Hamburg 

3^ now at Brown University, Providence, RI, USA 

3^^ supported by a MINERVA FeUowship 

now at ICEPP, Univ. of Tokyo, Tokyo, Japan 

present address: Tokyo Metropolitan College of Allied Medical Sciences, Tokyo 116, 
Japan 

^2 supported by the Pohsh State Committee for Scientific Research, grant No. 2P03B09308 



V 



supported by the Natural Sciences and Engineering Research Council of 
Canada (NSERC) 

supported by the FCAR of Quebec, Canada 

supported by the German Federal Ministry for Education and Science, 
Research and Technology (BMBF), under contract numbers 057BN19P, 
057FR19P, 057HH19P, 057HH29P 

supported by the MINERVA Gesellschaft fiir Forschung GmbH, the German 
Israeh Foundation, the U.S.-Israel Binational Science Foundation, and by the 
Israel Ministry of Science 

supported by the German-Israeli Foundation, the Israel Science Foundation, 
the U.S.-Israel Binational Science Foundation, and by the Israel Ministry of 
Science 

supported by the Italian National Institute for Nuclear Physics (INFN) 
supported by the Japanese Ministry of Education, Science and Culture (the 
Monbusho) and its grants for Scientific Research 

supported by the Korean Ministry of Education and Korea Science and Engi- 
neering Foundation 

supported by the Netherlands Foundation for Research on Matter (FOM) 
supported by the Pohsh State Committee for Scientific Research, grant 
No. 115/E-343/SPUB/P03/002/97, 2P03B10512, 2P03B10612, 2P03B14212, 

2P03B10412 

supported by the Polish State Committee for Scientific Research (grant No. 
2P03B08614) and Foundation for Pofish-German Collaboration 
partially supported by the German Federal Ministry for Education and Science, 
Research and Technology (BMBF) 

supported by the Fund for Fundamental Research of Russian Ministry for 
Science and Education and by the German Federal Ministry for Education 
and Science, Research and Technology (BMBF) 

supported by the Spanish Ministry of Education and Science through funds 
provided by CICYT 

supported by the Particle Physics and Astronomy Research Council 
supported by the US Department of Energy 
supported by the US National Science Foundation 



VI 



1 Introduction 



Supersymmetry (SUSY) theories relate bosons and fermions by associating to each fermion 
a bosonic partner and vice-versa. Among the appeahng consequences of this symmetry 
is the cancellation of quadratic divergences occurring in the scalar Higgs sector of the 
Standard Model (SM) and models beyond the SM 0. There is, on the other hand, no 
experimental evidence for SUSY. Since supersymmetric particles are not observed at the 
masses of their standard partners, SUSY would need to be broken. 

In the supersymmetric extension of the SM known as the Minimal Supersymmetric Stan- 
dard Model (MSSM), the supersymmetry-breaking terms are added by hand, generating 
a model with many free parameters (for a review and the original references see 0). In 
this model the breaking of SU(2) x U(1)y is generated through the vacuum expectation 
values f 1 and V2 of two Higgs doublets, which give masses to the down-type quarks and 
charged leptons {yi) and to the up-type quarks (^2). Selectrons (cl, cr) and squarks 
(^L, g^) are the scalar partners of the left- and right-handed electrons and the quarks 
of flavor f. The supersymmetric partners of the gauge bosons and the Higgs particles, 
known as gauginos and higgsinos, mix together giving rise to four neutral mass eigen- 
states Xi (neutralinos) and two charged mass eigenstates xf (charginos). The masses of 
the neutralinos depend on the supersymmetry-breaking parameters Mi and M2 (for the 
U(l) and SU(2) gauginos), on the higgsino mass parameter /i and on the ratio of the two 
Higgs vacuum expectation values tan/? = V2/V1. There is a new multiplicative quantum 
number, i?-parity {Rp), which takes the values +1 for the ordinary particles and -1 for 
the supersymmetric particles. In models where Rp is conserved, the supersymmetric par- 
ticles are produced only in pairs and the lightest supersymmetric particle (LSP) is stable. 
In these models, the production of a slepton and a squark is the lowest order process in 
which supersymmetric particles could be produced at HERA [Q, ||, ||. Other processes, 
such as the production of a slepton and a gaugino have a smaller cross section. 

We report here on the search for the i?p-conserving production and decay of a selectron Q 
and a squark, e~^p — *• e^qlX (a = L, R, b = L,R). This process is mediated by the t- 
channel exchange of a neutralino, as shown in figure 0. As the production of high-mass 
particles in the final state involves high Bjorken-x valence quarks from the proton, the 
process is mostly sensitive to the up quark. This measurement was performed with the 
ZEUS detector using an integrated luminosity C = 46.6 pb~^. The HI collaboration has 
published results of a search similar to the one presented in this paper, using an integrated 
luminosity C ~ 6.4 pb~^ 0]. 

The selectron (squark) can decay directly to the lightest neutralino, Xi^ ^iid a positron 
(quark): e ex? {<j —>■ qXi)- Under the assumption that the lightest neutralino is the 
LSP and that Rp is conserved, one can conclude that the neutralino will escape detection. 
In this case, the signature for the production of a selectron and a squark is one positron 
from e decay, a high Pt hadronic system from q decay and missing momentum from the 
two neutralinos [§. The search for this process is used to set limits on the masses of 

^Throughout this paper we wiU call "selectron" the scalar partner of the positron. 



the selectron and of the squark for a wide range of values of the MSSM parameters. 
Limits are derived for all squarks assuming that they are degenerate in mass. We also set 
limits separately for the u squark. These limits can be considered complementary to the 
strong limits on the squark mass from pp collider experiments |P that are valid within 
models with additional constraints compared to the MSSM. The present search leads to 
limits comparable to those so far obtained from e~^e~ experiments [10, 11 1 for particular 
combinations of the e and q masses. 



2 Experimental Setup 

The 46.6 pb^^ of data used for this analysis were collected with the ZEUS detector at 
the HERA collider during the years 1994 to 1997. The collider operated at a center-of- 
mass energy of 300 GeV with positrons of energy E^. = 27.5 GeV and protons of energy 
Ep = 820 GeV. 

The ZEUS detector is described in detail elsewhere ||T2|. The main subdetectors used in 
the present analysis are the central tracking detector (CTD) positioned in a 1.43 T 
solenoidal magnetic field and the compensating uranium-scintillator sampling calorimeter 
(GAL) g, subdivided in forward (FGAL), barrel (BGAL) and rear (RGAL) sections. 



Under test beam conditions the GAL energy resolution is 18%/ J E{GeV) for electrons 



and 35%/ y E{GeV) for hadrons. A three-level trigger was used to select events online. 
The trigger criteria applied relied primarily on the energies deposited in the calorimeter. 
The trigger decision was based on electromagnetic energy, total transverse energy and 
missing transverse momentum. Timing cuts were used to reject beam-gas interactions 
and cosmic rays. 

The luminosity was measured to a precision of 1.5% from the rate of energetic brems- 



strahlung photons produced in the process p ep7 [ITH]. 



3 Production Model 

The amplitude of the process depicted in figure |I| is given by the sum of four exchange 
graphs, one for each MSSM neutralino. The cross section depends on the MSSM pa- 
rameters Ml, M2, /i, tan/3, and on the masses of the produced particles, Cep^e^qlx = 
cr^b(Mi, M2, /i, tan/3, me^, m^f ). The cross section depends to a very good approximation 
on the sum mg + mq of the scalar particle masses. The branching ratios, B, for the decays 
e — > exi and q — >• gxi depend on the same MSSM parameters. Squarks are not allowed 
to decay to gluinos (g) as we assume rrig > rriq. 

For < Ml and M2, the x? is higgsino-like and the cross sections are small, due to the 
small coupling constants. In the region |/u| > Mi and M2, the Xi is gaugino-like and good 
limits on the sparticle masses are possible with the exception of particular combinations 
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of Ml and M2 (for > 0) when the LSP becomes a chargino. In the asymptotic region 
\fi\ ^ M2 > Ml, the mass of the hghtest neutrahno is m^o ~ Mi and the masses of 
the next-to-hghtest neutrahno and of the charginos are m^o ~ m + ~ M2; therefore for 
M2 > mg and rriq, the scalar particles can only decay to the LSP and B = I, while 
B ^ 30% when M2 < mg and rriq. In this region the cross section depends weakly on 
tan/3. We take = 500 GeV as an example of large /i values. 

To reduce the number of free parameters, the cl and cr are assumed to have the same mass 
mg, and all the squarks (except i whose contribution can be neglected) to have the same 
mass rriq. Alternatively we consider the case that u is the only squark contributing to the 
process. We assume no mixing between cl (^l) and cr (g^). In Grand Unified Theories 
(GUT) the gaugino masses unify at the GUT scale which leads, at the electroweak scale, 
to the relation Mi = |tan^6'vyM2 ~ |^2, which is adopted here to fix the branching 
ratios and compare with other experiments. 



4 Monte Carlo Simulation 



Monte Carlo simulations are used to determine the efficiency for selecting the signal, 
and also to estimate the rate of the SM background. All the signal and background 
events were passed through the simulation of the detector response based on GEANT 
1^], incorporating the effects of the trigger. They were subsequently processed with the 



same reconstruction and analysis programs used for the data. 



4.1 Signal 



The calculation of the leading-order cross sections and the simulation of the signal have 
been performed with a Monte Carlo (MC) generator (MSSM, see [|17|) based on the cal- 



culations of 0. In the calculation of the total cross section, the proton parton density is 
modeled using the GRV-LO [T^ set. The change of the structure functions to CTEQ-LO 
[p!9|] and GRV-HO produced variations always less than 4%. The effect of the initial state 
QED radiation from the incoming lepton was included using the Weizsacker- Williams 
approximation. To generate the signal MC samples the MSSM MC is interfaced to ARI- 
ADNE [|0l, which implements the Color Dipole Model (CDM) to simulate the QCD 



radiation, and JETSET |]2T]], for the fragmentation of the final state partons. The signal 
samples were generated for 118 different combinations of the values for the masses mg, rriq 
and m^o, by assuming the exchanged neutrahno to be a pure photino (7), since the event 
kinematics and the distributions of the decay products are insensitive to the neutrahno 
mixing parameters. The e and q were then forced to decay directly into Xi- 
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4.2 Background 



To evaluate the SM background we considered the following reactions which can mimic 
the signature for a selectron-squark event: 



• neutral current (NC) deep inelastic scattering (DIS) events with missing momen- 
tum due to the presence of neutrinos, muons or not-completely-contained hadronic 
energy; 

• charged current (CC) DIS events with a true or fake positron in the final state; 

• W production processes which can mimic the signal via the leptonic decays; 

• lepton pair production from Bethe-Heitler interactions (77 e^e~ , ,t~^t~); 

• photoproduction with high transverse energy and heavy flavor production, with true 
or fake positrons in the flnal state. 



Backgrounds from NC and CC DIS events were simulated using HERACLES ||2^, which 
includes first order electroweak radiative corrections. The hadronic final states were sim- 



ulated using LEPTO [|3[ and JETSET interfaced with HERACLES using DJANGO p4 
QCD radiation was simulated using the CDM as implemented in ARIADNE and 



as 



a cross check, with the exact first order matrix elements followed by parton showering 
(MEPS option). The production of was simulated using EPVEC ^E\. The produc- 



tion of lepton pairs in both elastic and inelastic Bethe-Heitler interactions is simulated 
using LPAIR ||2^. All these samples correspond to luminosities at least five times larger 
than that of the data. The photoproduction background was simulated using events with 
high transverse energy produced in direct- and resolved-photon interactions, which were 
generated using HERWIG [^]. Samples of heavy fiavors (cc, bb) were also produced via 
the boson gluon fusion mechanism, as implemented in AROMA pS) . These last two sam- 
ples, which represent a small background, correspond to a luminosity comparable to that 
of the data. 



5 Event Selection 



The event selection for the SUSY particles search requires the following global event 
quantities calculated from the reconstructed vertex and the calorimeter measurements: 

i 
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where E[ is the energy measured in the i^^ calorimeter cell and Pi = E^fli, with ni the 
unit vector from the reconstructed vertex to the cell center^. The sums run over all 
calorimeter cells. The missing transverse momentum ft and the longitudinal momentum 
variable E — Pz equal, respectively, zero and twice the incident positron beam energy 
{2Ee = 55 GeV) for fully contained events. The signal events are expected to have large 
/t and low E — Pz due to the momentum carried away by the escaping neutralinos. The 
calculation of the hadronic transverse momentum (P^) excludes cells in the calorimeter 
sum which belong to a cluster which is identified as the candidate positron from e decay. 

Events are selected by requiring an identified positron with transverse momentum P^ > 4 
GeV and a hadronic system with P^^ > 4 GeV. 



Positron candidates were identified using the pattern of energy deposits in GAL |29|, which 
provided the measurement of its energy Ee and angles 6e, (pe- Additional requirements 
were applied: 



• isolation: the total transverse energy not associated with the positron candidate 
within a cone of radius R=0.8 units in the {t], 0) plane around its direction must 
be less than 3 GeV and the sum of the momenta of the tracks pointing in a cone 
of radius R=0.8 units in ri,(j) around the positron, excluding the highest momentum 
track, must be less than 2 GeV. These cuts remove fake positrons from GG and 
photoproduction that are close to the hadronic jet; 

• track matching: if the candidate positron is within a region where the tracking 
efficiency is high, namely 6e > 0.35 radians, a match with a track is required such 
that 

— the distance of closest approach between the electromagnetic (EM) cluster and 
the track extrapolated from the GTD be less than 6 cm, 

— the ratio of the track to GAL transverse momentum, p^t^ack^pe ^ g 2^ 

— the track has a positive charge or p^'^'^'^^ > 45 GeV; 

if 6e < 0.35, the requirement on the candidate positron transverse momentum is 
raised to P^P > 10 GeV; 

• fiducial cut: polar angle 6e < 2; this cut eliminates most of the low NG DIS 
background. 



To ensure that the selected events come from e~^p interactions we apply standard cosmic 
ray and beam gas rejection cuts and require a reconstructed vertex within 50 cm of the 
nominal interaction point along the Z axis. 

^The ZEUS coordinate system is a right handed system with the Z axis pointing in the proton direction. 
The pseudorapidity is defined as = — ln(tan |), and the polar angle 9 is measured with respect to the 
Z direction. 
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In figures | a) and b) the ft and E — Pz distributions of the data, after a cut ft > 10 GeV, 
(full points) are compared with the total background expectation (full line) and with a 
Monte Carlo signal sample corresponding to the masses mg = 80 GeV, rrig = 80 GeV 
and m^o = 50 GeV (dotted line). Good agreement is observed between the distributions 
of the data and the background Monte Carlo. 

To determine the values of the final cuts, an optimization procedure was applied, which 
maximized the ratio e/U for simulated events, where e is the efficiency and U is the 95% 
confidence level (CL) upper limit on the signal when the number of events found equals the 
expected background. Events were accepted if they satisfied the following requirements: 

• E-Pz< 50 GeV; 

• ft> 14 GeV; 

. (E-Pz)/ft< 1. 

The acceptance for the e+ q Xi Xi final state is close to zero for small mass differences 
Am = min(me — m^o, m^ — m^o) and reaches a plateau for Am > 10 GeV. The level of the 
plateau increases from 25% at rrig = rrig = 40 GeV to about 50% at m^ or rriq = 120 GeV. 
The acceptance was modeled with a two parameter fit as a function of m^o for fixed m^ 
and rriq. Subsequently a bilinear interpolation was used to obtain the parameters at any 
mg and rriq. 

In figures c), d) and e), the data, the signal sample and the SM background are shown 
in the ft,E — Pz plane, where the final cuts are also drawn. The expected SM back- 
ground is I.QQIqII events, as shown in table [|. The systematic uncertainties quoted for 
the background are evaluated as described in section |^. No photoproduction or heavy 
flavor event survived the selection. An alternative method was used to confirm that the 
photoproduction background was not underestimated due to the limited MC statistics. 
It was found that the rate of photoproduction events passing all of the kinematic cuts 
(except positron identification) was only a small fraction of the CC background rate. As- 
suming that the fake positron fraction was the same for CC and photoproduction events, 
the resulting photoproduction rate was negligible. 

One event survived the selection criteria and was identified as containing a high-Q^ 
positron with associated f in the calorimeter due to two muons in the final state. 

6 Systematic Uncertainties 

The main source of systematic uncertainty on the efficiency e comes from the interpolation 
between the generated Monte Carlo mass points (Ae/e = 6%). Other contributions from 
detector effects that affect both the acceptance and the background are the uncertainties 
on the positron identification efficiency (Ae/e = 1%), on the calorimeter energy scales 
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(Ae/e = 2%), on the efficiency of the matching between a track and an EM cluster 
(Ae/e = 3.5%), on the isolation requirements (Ae/e = 3%) and on the vertex distribution 
(Ae/e = 1%). Another common source is the 1.5% uncertainty in the total integrated 
luminosity of the 1994-1997 e^p data sample. We quote an overall systematic uncertainty 
Ae/e = 8.5%. 

In addition to the experimental sources, effects due to the different treatment of the QCD 
radiation in the final state for the CC and NC background were evaluated using Monte 
Carlo samples generated with the MEPS option. The variation of the CC background 
was found to be below the statistical uncertainty. No NC MEPS event survived the final 
selection. This difference with respect to the ARIADNE NC background is included as 
a systematic uncertainty. The theoretical uncertainties on the cross sections of the back- 
ground processes are also included in table |I[ The quoted errors are due to uncertain- 
ties on: the proton parton distribution functions (pdf) for the NC and CC backgrounds 
(5 — 6%), the missing higher orders and proton and photon pdf for W production (35%) 
and the uncertainties mainly due to inelastic contributions on the lepton pair production 
(10%). The total error on the background (~ ±35%) is taken to be the sum in quadrature 
of the statistical and systematic errors. For the estimated background and one observed 
event this ~ 35% variation of the background normalization results in a 4% variation of 
the 95% CL upper limit on the number of signal events. 

7 Results 

The event found in the data sample is compatible with the expected SM background. 
Therefore we derive upper limits on the e,q production cross section times the branching 
ratios [a x B) for the decay to the lightest neutralino Xv The 95% CL upper limit on 
the signal, including the systematic uncertainty on the acceptance and the uncertainty in 
the background as described in ||3^, is A^ul = 3.9 events. The 95% CL upper limit on 
cr X 5 is obtained as a function of the selectron, squark and neutralino masses through 
a X B < N\ji^/eC, where e = e{me,mq,m^o) is the parametrized efficiency and C is 
the integrated luminosity. Figures ^ a) and b) show the upper limit on cr x i? in the 
plane defined by mg and rriq for two values of m^o (35 and 50 GeV). At high mg and 
niq and large Am we exclude a x B > 0.21 pb. For m^o < 50 GeV, Am > 10 GeV and 
(mg + mq)/2 > 70 GeV we exclude a x B > 0.25 pb. For any mass combination with 
Am > 10 GeV we exclude a x B > 0.5 pb. 

From the comparison of the upper limit on a x B with the theoretical value from the 
model calculations we derive exclusion areas in the parameter space of the MSSM. The 
theoretical cross sections were evaluated by taking into account the complete neutralino 
mixing and the branching ratios for the decays of e and g to Xi? computed at the tree 
level. The cross section times branching ratios is computed from: 

(a X B\^ = E E { ^e+p^..,^ 5(e+ - e+x?) B{qi ^ g^x?) + 

a,b=L,R f=l 
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In figures ^ c) and d) we show the excluded regions in the plane defined by the lightest 
neutralino and by (mg + mq)/2. The cross section for the e and q production depends 
mainly on the sum of the two masses; we give the excluded region for mg = mg since the 
efficiency and the branching ratios depend separately on the two masses. The limits are 
shown for /x = ±500 GeV as examples of large values of /x. The limits for = +500 GeV 
and /i = —500 GeV are similar and differ by the hatched region in the figures. Limits 
are also shown for the intermediate value /i = —100 GeV. For large the excluded 
region reaches (mg + mq)/2 = 77 GeV for a 40 GeV neutralino. This limit worsens at 
lower neutralino masses, because new decay channels to charginos and next to lightest 
neutralino open and compete with the direct decay to Xv the limit M2 ^ Mi the 
charginos and the next to lightest neutralino masses increase leaving only the direct 
decay channel to Xi open. In this case the excluded region reaches (mg + mq)/2 = 84 
GeV for a massless neutralino. The excluded region is limited by the small cross section 
of the process at large (mg + rrig) /2, while for large neutralino masses it is limited by the 
efficiency that falls to zero as Am — > 0. A large variation in tan/5 results only in slight 
changes. The limits are shown in c) for tan/? = 1.41 and in d) for tan/? = 10 as examples 
of a small and a large value of tan/?. The previous limits obtained by HI are also shown 
in figure ^ c). 

Figure ^ shows the excluded regions in the plane rrie, rriq for fixed values of m^o = 35, 50 
GeV and for different combinations of n and tan (3. The excluded area is approximately 
triangular, defined by a line of constant mg+mg and by two lines parallel to the coordinate 
axes due to the low efficiency for Am < 10 GeV. In figures § a) and b) the limits are 
given for the large /z region (|/i| = 500 GeV) and tan/5 =1.41, 10. The limits obtained at 
tan/3 = 1.41 and fi = —100 GeV are shown in figure ^ c). 

In the range 45 < (mg + mg)/2 < 85 GeV the up quark contribution to the cross section 
ranges between 70% and 90% because it dominates the parton densities at high-x. The 
limits on the u mass, assuming all the other squarks to be much heavier, are only ~ 2 
GeV below the limit obtained for degenerate squark masses. This is shown as an example 
in figure H d), to be compared with figure § c). 



8 Comparison with other experiments 

A search for the production of selectrons and squarks at HERA has been published by HI 
[0. The present analysis substantially improves these limits due to the seven-fold increase 
in the integrated luminosity. 

Published limits from LEP at ^/s = 161 — 172 GeV exclude selectrons with masses lower 



than ~ 80 GeV |jTO[. The present analysis goes beyond this limit for squark masses below 
~ 80 GeV. LEP experiments have also reported limits on the stop and on the sbottom 
squarks [|ll|] of m^, m^ ~ 75 GeV. The limits given in our analysis are (for mg ~ rrig) at 



the same level as those obtained by combining the LEP limits on mg with those on m 
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and mr extended to the other squarks, including u. Present hmits on the mass of the 



hghtest neutrahno |^ exclude m^o < 25 GeV for tan /? ~ 1 and m^o < 40 GeV for large 
tan p. 

Strong limits on the squark mass have been obtained from the pp experiments at the 
Tevatron and are complementary to those obtained in the present analysis. The 
Tevatron experiments adopt strong cuts on /t and select final states with many high- 
Pt jets or leptons to exclude the QCD background. Therefore, they are more sensitive 
to squark cascade decays than to the direct decay q — >■ qxi, especially when the mass 
difference rrig — m^o is small (~ 20 GeV). On the other hand, the search described in this 
paper looks only for direct decays and, due to the relatively low cuts on ft (14 GeV) and 
Pf^ (4 GeV), is sensitive to smaller mass differences. Limits on SUSY are often expressed 
in terms of the minimal supergravity model (mSUGRA) in which the masses of scalars 
and fermions are no longer independent parameters [Q. In this model the squarks are 
always heavier than selectrons and gluinos, and the dominant squark decay is to gluino and 
quark. In the region where the squark is lighter than the gluino, the Tevatron experiments 
give limits under more general assumptions than mSUGRA but keep the GUT inspired 
assumption that links the masses of the neutralinos and the charginos to that of the 
gluinos. In this framework limits cannot be set on light squarks (~ 80 GeV) when the 
gluino is too heavy (~ 500 GeV), due to the lack of sensitivity at small nig — m^o. 

The present analysis investigates a different region of the MSSM parameter space, where 
the gluinos are heavier than the squarks and the squarks masses are independent of the 
selectron masses, allowing also for rriq < mg- Moreover, no relationships are assumed 
between m^o and rrig. 



9 Conclusions 



We have searched for the SUSY process e^p — >• e~^qX (e"*" — e^Xi ; q ~^ qXi) using 
46.6 pb~^ of data at an e^p center-of-mass energy of 300 GeV recorded with the ZEUS 
detector at HERA. One candidate event is found while 1.99to;84 are expected from Stan- 
dard Model processes. The upper limit on the production cross section times branching 
ratios is found to be cr x i? < 0.5 pb at the 95% CL for mass differences mg — m^o and 
rriq — m^o greater than 10 GeV. Excluded regions in the MSSM parameter space have 
been derived. We exclude (mg + mq)/2 < 77 GeV at the 95% CL for m^o = 40 GeV and 
large values of the MSSM parameter |yu|. The process is dominated by the u contribution 
and the exclusion limit is 75 GeV when only the u squark is considered. 
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Process 


Expected Events 


Stat. 


Theor. 


Syst. 


Total 


NC DIS 


0.60 


±0.29 


±0.03 


+U.U6 

-0.67 


+U.3U 
-0.73 


CC DIS 


0.18 


±0.07 


±0.01 


+U.24 
-0.10 


+0.25 
-0.12 


W 


1.04 


±0.04 


±0.36 


+u.uy 

-0.12 


+0.37 
-0.38 


/+/- 


0.17 


±0.12 


±0.02 


+U.17 
-0.08 


+0.21 
-0.14 


Total 


1.99 


±0.32 


±0.36 


+U.3i 

-0.69 


+0.57 
-0.84 



Table 1: Background expectation from MC calculations. The table shows the background 
process, the expected number of events normalized to data luminosity, the statistical 
error related to the MC statistics, the error due to the theoretical uncertainty on the 
cross section, the systematic error and the total error. 




Figure 1: Diagram for the selectron and squark production process. 
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Figure 2: Distributions of events with a selected positron and /t> 10 GeV, in a) and 
b) E — Pz for data (points), SM background (full line) and a signal example (mg = 80 
GeV, nig = 80 GcV, m^o = 50 GeV) (dashed line). The distributions of events in the 
E — Pz versus /t plane are shown in c) for the data, in d) for the MSSM example and in e) 
for the SM background, where the MC samples are normahzed to 5 times the luminosity 
of the data. 
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Figure 3: a) and b) 95% CL upper limits on cr x S in the plane defined by nig and mg 
for two values of the mass of the lightest neutralino: m^o = 35 GeV in a) and m^o = 50 
GeV in b). c) and d) show the regions excluded at the 95 % CL in the plane defined by 
(mg + mq)/2 and Tn^o, where limits are evaluated along mg = niq for tan/3 = 1.41 (c)) 
and for tan/3 = 10 (d)). The grey area is for = 500 GeV. For /i < the excluded 
region covers also the hatched grey area. The dashed line is for /i = —100 GeV. These 
limits are obtained for M2 = 2Mi. If this relation is modified to M2 ^ Mi the excluded 
area for \n\ = 500 reaches the full line. The kinematic limit is indicated by the straight 
line in the upper-left corner. Previous HI limits from [|^ are also shown in c). 
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Figure 4: Excluded regions at the 95% CL in the plane defined by the selectron and 
squark mass, for m^o = 35 GeV (grey area) and m^o = 50 GeV (double-hatched area). 
In a) (tan/3 = 1.41) and b) (tan/? = 10) the limits are for = 500 GeV. For < the 
excluded region includes also the single-hatched area. The limits obtained for fj, = —100 
GeV and tan = 1.41 are shown in c), where LEP limits on mg are also given. The limits 
for the up squark alone, for the same values of /i and tan/5, are shown in d). 



15 



